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ABSTRACT 

The propagation of cosmic -ray protons in the Galaxy is discussed under the framework of 
a three dimensional convection-diffusion model. Starting with the assumption of a uniform 
and continuous distribution of cosmic-ray sources injecting CRs continuously in the Galaxy 
and by invoking a supernova explosion at various distances from the Earth, it is found that 
only those sources located within a distance of ~ l.bkpc can produce appreciable temporal 
fluctuations in the CR proton flux observed at the Earth. So, the construction of the local CR 
proton spectrum is discussed by seperating the contributions of the distant sources from that 
of the nearby sources. The contribution from the distant sources is treated in the framework 
of a continuous source distribution model both in space as well as time, but that of the nearby 
sources in a discrete space-time source model. The study predicts the presence of at least one 
old nearby source with a characteristic age of ~ 10 5 yrs located at a distance of ~ Q.lkpc to 
explain the observed proton flux below ~ lQOGeV. 
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1 INTRODUCTION 



Considering the similarity of the power supplied by a super- 
nova explosion (SNE) which is approximately 10 51 erg /30yrs ~ 
10 42 erg/ s to the power required to maintain the cosmic-ray (CR) 
energy density in the Galaxy (~ 10 41 erg/ s), it is widely beleived 
that the majority of CRs upto ~ WOTeV are accelerated in super- 
nova remnant (SNR) shock waves (Blandford & Eichler 1987). At 
least for CR electrons, this hypothesis is strongly supported by the 
detections of nonthermal X-ray synchrotron emission from SNRs 
like SN 1006 (Koyama et al.1995), RX J1713.7-3946 (Koyama et 
at. 1997), Cas A (Allen et at. 1997), IC 443 (Keohane et at. 1997), 
G374.3-0.5 (Slane et al. 1999) and RX J0852.0-4622 (Slane et al. 
2001). But, for CR protons, their acceleration in SNRs is still not 
fully understood (Berezhko et al.2002). 

Several models have been proposed to discuss the CR propa- 
gation in the Galaxy. Among these models, for CR electrons the 
idea of separation of distant and nearby sources can be found 
in Shen (1970), Shen and Mao (1971), Cowsik and Lee (1979), 
Atoyan et al. 1995 and Kobayashi et al. 2004. For CR protons and 
nuclei, this approach had been discussed by Lingenfelter (1969) 
using a simple model considering only the particle losses due to 
escape. 

Much later, Erlykin & Wolfendale (2001a, 2001b) studied the 
CR spectra at the Earth using Monte Carlo calculations assuming 
SNRs as the main sources of CRs. In their model, CR nuclei are 
assumed to be accelerated throughout the SNR lifetime and liber- 
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ated into the ISM when the acceleration process ceases. They as- 
sumed the SN to be distributed stochastically throughout the nearby 
Galaxy. Their CR propagation model in the Galaxy includes only 
the diffusive process, neglecting convection due to the Galactic 
wind, energy and catastrophic losses due to interactions with the 
ISM. 

Strong & Moskalenko (2001) have studied the influence of the 
discrete nature of SNRs on the CR protons spectra in the Galaxy us- 
ing their galprop CR propagation code. Their calculations involved 
the stochastic nature of SNR events in the Galaxy and a realistic 
gas distribution for CR interactions. 

In this paper, the contributions of the distant and nearby 
sources to the local CR proton flux are treated seperately. The 
study considers only the observed local SNRs as the nearby sources 
unlike in earlier works and the propagation of CR protons in 
the Galaxy is treated under the framework of a 3-D convection- 
diffusion model in which the protons undergo convection, diffu- 
sion, adiabatic deceleration, ionization and coulomb losses, parti- 
cle loss due to pion production processes and injection from the 
sources. The calculations are done mainly for the CR flux in the 
vicinity of the Sun. For the separate treatment of the distant and 
nearby sources, CRs from the distant sources are treated in the 
framework of a continuous source distribution model both in space 
and time while those from the nearby sources are considered in a 
discrete (in space and time) source model. Although the model pre- 
sented here is true for any type of pointlike sources, the calculations 
are done considering SNRs as the CR sources. 

The study given in this paper starts with the assumption that 
the CR sources distribution in the Galaxy are spatially uniform and 
temporally continuous and then, by invoking a supernova explo- 
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sion at various distances from the Earth, it is found that only those 
sources located close to the Earth can produce appreciable temporal 
fluctuations in the local CR proton flux. Knowing this and consid- 
ering only the observed local SNRs, the study predicts the presence 
of atleast one old nearby source (which is referred to as the missing 
source in this paper) in order to explain the observed proton flux 
below ~ WQGeV. 

The plan of this paper is as follows. Section 2 discusses about 
the propagation of CR protons in the Galaxy as well as the spatial 
and temporal dependence of their flux from a single source. Section 

3 studies their number density variations in the Galaxy and section 

4 discusses about the CR proton flux in the Solar vicinity. Finally, 
section 5 presents the results and discussions about the study pre- 
sented in this paper. 



2 PROPAGATION OF CR PROTONS IN THE GALAXY 
2.1 Propagation properties of CR protons 

CRs undergo diffusion in the Galaxy due to scattering either by 
magnetic field irregularities or by self excited Alfven and hydro- 
magnetic waves. The diffusion coefficient in the Galaxy is gener- 
ally obtained using the observed secondary to primary B/C ratios. 
For energies upto around lOOGeV^, it has the form D(E) oc E 
with S = 0(0.6) for E < 5GeV(^ 5GeV) (Kobayashi et al. 
2004 and references therein). However, in the TeV region there is 
a strong argument from the studies of CR anisotropy that a value 
of S ~ 0.6 gives an anisotropy value much larger than the ob- 
served value of ~ 10~ 3 (Ptuskin 2005, Ambrosio et al. 2003 and 
references therein). Also, experimental studies on B/C and (sec- 
ondary species of Fe)/Fe had supported the idea that S changes 
from 0.6 - 0.3 in the GeV - TeV region (Furukawa et al. 2003). 
But, these constraints imposed by the observed anisotropy are not 
fully well defined because of the presence of solar wind, magnetic 
field discontinuities in the solar cavity, large scale structure of local 
magnetic field lines and the measurements of only the CR velocity 
component perpendicular to the rotation axis of the Earth (Cesarsky 
1980). There is also a difficulty related to the proton spectrum in 
changing 8 from 0.6 to 0.3 since the proton spectrum would then 
have a break at around ITeV for a power law injection spectrum 
whereas observations indicate that the proton spectrum follow a 
continuous single power law throughout the (10 — 10°)GeV en- 
ergy range. Therefore, in this paper the diffusion coefficient is as- 
sumed to follow a single power law in the high energy region and 
takes the form (Atoyan et al. 1995) 



D(E) = D 1 + 



E 
E~ 



(1) 



where Eq = 3GeV, a — 0.6 and Do is chosen such that 
D(WGeV) = 10 28 cm 2 s~\ Eq. (1) is considered to be valid 
atleast upto ~ lOOTeV which is approximately the maximum en- 
ergy of accelerated particles by SNR shock waves (Lagage &: Ce- 
sarsky 1983). 

While propagating in the Galaxy although high energy CR 
protons do not suffer from radiative energy losses like the high en- 
ergy CR electrons, protons with energies less than around lGeV 
suffer significant energy losses due to ionization and Coulomb in- 
teractions with the interstellar medium (ISM). Their energy loss 



rate is given by (Mannheim and Schlickeiser 1994) 
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where E is in GeV, m p — 0.938GeV is the proton mass en- 
ergy, uh = Uhi + 2nH 2 is the ISM hydrogen atom density 
and linn is the density of ionised gas in the ISM (Cordes et al. 
1991). The atomic hydrogen density uhi is taken from (Gordon 
& Burton 1976) and the molecular hydrogen density uh 2 from 
(Bronfman et al. 1988). For the present study, their values are cho- 
sen at the position of our Solar system as nu = l.llcm~ 3 and 
Uhii — 0.022cm -3 . In propagation models like the one which 
is going to be presented here which include a Galactic wind in the 
direction perpendicular to the Galactic plane, protons with energy 
greater than ~ lGeV lost their energy mainly due to adiabatic 
cooling if the wind velocity depends on the Galactic height. In 
addition to this, CRs are also convected away from the Galactic 
plane due to the Galactic wind. Thus low energy CRs with energy 
less than ~ 100 GeV are not able to reach the Earth effectively 
from distant sources because of their slow diffusion since diffusive 
timescale goes as tdiff oc 1/D(E). Hence, the effects of vari- 
ous energy loss processes and the Galactic wind are important for 
low energy protons and should be considered while discussing their 
propagation in the Galaxy. 



2.2 CRs from a single point-like source 

In the framework of a convection- diffusion model, the propaga- 
tion of cosmic ray protons in the Galaxy can be represented by the 
equation 



V • (DX7N - \N) + A 



V- V 



J 'pp 
„ dN 



(3) 



where E is the proton kinetic energy in GeV, N(r, E) is the dif- 
ferential number density , D(r,E) = D(E) is the diffusion co- 
efficient which is assumed to be spatially uniform, Q(r,E,t) is 
the proton production rate, i.e. Q(r, E,t)d 3 rdEdt is the num- 
ber of protons produced by the source in a volume element 
d 3 r = dxdydz in the energy range (E, E + dE) in time dt, 
V(r) = SVozz is the convection velocity in the direction per- 
pendicular to the Galactic plane with Vo = 15 km s~ kpc" 1 
and b(r, E) = —B(E) is the proton energy loss rate due to ion- 
ization and Coulomb interactions which is assumed to be spa- 
tially constant. [—N p (E,r)/T pp (E,r)], r pp — E/(dE/dt) w and 
(dE/dt)^ — 8 x 10^ 16 hhE in GeV s _1 are the catastrophic loss 
term, energy loss time scale and energy loss rate respectively due 
to pion production processes. With these substitutions, Eq. (3) can 
be simplified in rectangular coordinates as 
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The Green's function G(x, x ,y,y ,z,z ,E,E ,t,t ) of Eq. (4), 
i.e. the solution for a ^-function type source term Q(x,y, z, E, t) = 
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5(x — x )8(y — y )S(z — z )8(E — E )8(t — t ) is obtained as 
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where P{E) = V E + B(E), A pp = 8 x 10 _16 n H and 

du 



A(E) — exp 



Wo 



V u + B(u) 
Then, the general solution of Eq. (4) is given by 



N(x,y,z,E,t) = 



dx 



dy 



dz / dE 

IE 



dt Q(x ,y ,z , E ,t ) 

xG(x,x',y,y,z,z,E,E,t,t') (6) 

Now, to calculate the number density N(x,y,z,E,t) due to a point 
source located at (0, 0, 0) producing q(E) protons of energy E per 
unit energy interval at the time of explosion, let us start by consider- 
ing a small source volume p bounded by \x\ — xi, \y\ = j/i, \z\ — 
z\ where protons are generated. Before the explosion protons are 
assumed to be distributed uniformly in the source volume with den- 
sity qo(E). In other words we can write 



inside the source, 
outside the source, (7) 



Q(x,y,z,E,t) = qp(E)q(t), 

= 0, 

and for a point source 
]imq (E)p = q(E) 



With this assumption the proton density outside the point 
source positioned at (0, 0, 0) can be obtained by introducing Eqs. 
(5)&(7) in (6) as 
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Assuming burstlike injection of CR protons from the source ex- 
ploding at time t = to, we can write q(t ) = 8(t —to) which 
gives 



N(x,y,z,E,t) = 
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In cylindrical co-ordinates (r,6, z), Eq. (9) becomes 
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It should be mentioned here that the proton source spectrum q(E) 
is taken as 



q(E) = k(E 2 + 2Em p ) 



-(r+i)/2 



(E + m P 



(11) 



The constant k is obtained by normalizing the proton source spec- 
trum such that the total amount of kinetic energy contained in the 
injected protons is 0.2 times the total supernova explosion energy. 

Fig. 1 shows the E 2 times the proton flux 1(E) = 
(c/4n)N(E) at a distance r = O.lfcpc, z — from a source at 
different times t after their burstlike injection into the ISM at time 
t — 0. It can be seen that the energy spectra are modified consid- 
erably at different times due to the combined effects of different 
propagation and energy loss rates of protons at different energies. 
At shorter times t < W 4 yrs, only the high energy particles are 
able to travel the distance because of their faster diffusion in the 
Galaxy. Low energy particles are not able to cover the distance in 
the short time. But at longer times t > W 4 yrs, the low energy 
particles also start reaching the observer and by that time the high 
energy particles has already passed by and hence, their flux is sig- 
nificant suppressed. At still longer times t > W 6 yrs, even the low 
energy protons have already passed by and the particle fluxes are 
found to decrease with time at all the energies. 

The particle fluxes are also found to depend strongly on the 
source distance r. The distance dependence is found to follow 
a Gaussian function (see Eq. 10). Since the CR sources are dis- 
tributed within a region very closed to the Galactic plane (half 



4 Thoudam Satyendra 



— " 0.001 




Kinetic energy E (GeV) 

Figure 1. Temporal behaviour of CR protons at a distance r = 
O.lkpc, z = from a pointlike source. The spectra numbered as 
(1, 2, 3, 4, 5, 6, 7, 8, 9) are calculated at times t = (10 3 , 3 X 10 3 , 10 4 , 3 X 
10 4 , 10 5 , 3 X 10 5 , 10 6 , 3 X 10 6 , 10 7 )yrs respectively. The calculations 



are done for total output source energy W 
0.6,D(10GeV) = 10 28 cm 2 s- 1 and V -- 
a burstlike injection of the protons at time t - 



10 bl ergs,T 



lbkms 
0. 



^kpc 



2.4, a = 
assuming 




Kinetic energy E (GeV) 

Figure 2. Various CR proton flux at the position of the Sun due to a single 
source located at various distances r. The spectra marked as (1, 2, 3, 4, 5) 
are calculated for source distances r = (0.05, 0.1,0.5, 1.0, 1.5)kpc at 
their respective times t m (E) = r 2 /6D(E) for proton energy E = 
10GeV, where t m is defined as the time at which the maximum flux of 
CR protons of energy E is observed after their injection from the sources 
at time t = 0. The model parameters are the same as in Fig.l. The Solar 
modulation parameter <E> = 400MV. 



thickness ~ 150pc) and also our Solar system is only 15pc away 
from the plane (Cohen 1995), it will be assumed from now on 
that the sources as well as the Sun are positioned at z = 0. Fig. 
2 shows the various CR proton flux in the solar vicinity due to 
a source located at various distances from the Sun. The maxi- 
mum flux of protons of energy E from a source after their injec- 
tion into the ISM at t = would be observed at time t m (E) = 
r 2 /6D(E). The fluxes shown in the figure are calculated for r — 
(0.05, 0.1, 0.5, 1.0, 1.5)fcpc at their respective i m 's for proton en- 
ergy E = WGeV. The magnitude of the maximum flux at WGeV 
is found to decrease with the source distance r. The calculation 
includes the effects of Solar modulation using the force field ap- 
proximation (Gleeson & Axford 1968) with modulation parameter 
$ = 4O0MV. 



2.3 CRs from the source distribution in the Galactic disk 

The propagation of CR protons in the Galaxy in the case of station- 
ary injection from a continuous and uniform source distribution in 
the Galactic disk can be represented by the equation 



V • (DViV - VAT) ■ 



_d_ 

BE 



V- V. 



E-b) N\ - 



N 



+Q = (12) 



where the symbols have the same representation as in section 2.2. 
The proton density at a given location in the Galactic disk due to all 
the sources located beyond a distance r from the point is obtained 
by solving Eq. (12) assuming a uniform and continuous distribution 
of sources in the disk and is given by (Thoudam 2006) 
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Here, the source spectrum q(E) is taken as the same form as in 
Eq. (11) with normalization k (in place of k) which is obtained 
assuming a Galactic (SNE) rate of 25 SNE Myr~ 1 kpc~ 2 in the 
disk (Grenier 2000) and 0.2 times the total SNE energy going into 
the proton kinetic energy. It should be noted that in the case of a 
stationary injection of CRs into the Galaxy from sources continu- 
ously distributed in the Galactic disk, the equilibrium CR proton 
spectrum in the disk depends weakly on the actual geometry of the 
source distribution for a constant CR energy density in the Galaxy. 



3 CR PROTON DENSITY VARIATIONS IN THE GALAXY 

Let us try to examine the sources which can produce appreciable 
temporal density fluctuations of CR protons at a given location in 
the Galaxy. The study starts with the assumption of a uniform and 
continuous distribution of cosmic-ray sources injecting CRs contin- 
uously in the Galaxy. These sources will contribute the background 
CRs in the study. Then, let us try to invoke a SNE at increasing 
distances from the point and see its effects on the CR density at the 
point. 

Figs. 3&4 shows the temporal variation of the proton den- 
sity at lOGeV and WTeV energies respectively for a SNE at 
r — (0.1, 0.5, 1, 1.5)kpc distances from the point. The calculation 
assumes a particle release time of t = from the source. It can be 
seen that the temporal density variation decreases as the distance of 
the SNR increases. For WGeV protons, the amplitude of variation 
decreases from ~ 90% at r = O.lkpc to ~ 0.01% at r = 1.5kpc. 
Similarly for WTeV protons, the amplitude varies as ~ 700% to 
~ 0.2% as d increases from O.lkpc to 1.5kpc. The higher am- 
plitude in the case of lOTeV protons is expected because of the 
energy dependence of the diffusion coefficient which allow higher 
energy particles to diffuse faster than the low energy ones in the 
Galaxy. With this study, it can be concluded that the propagation of 
CRs from nearby sources (say within ~ 1.5fcpc) should be treated 
in a time dependent model whereas those from far sources (those 
beyond a distance of ~ 1.5fcpc) can be discussed in a steady state 
model. Also, the spatial treatment between the nearby and distant 
sources should be different since the mean distance between the 
sources can be considered to be comparable to their distances from 
the Earth for nearby sources whereas that of the distant sources can 
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Table 1. Parameters of SNRs with known ages located within a distance of 
l.hkpc from the Sun : 



SNR name 


distance r 


Age* 




(kpc) 


(yrs) 


G65.3+5.7 


1.0 


14000 


G73.9+0.9 


1.3 


10000 


Cygnus Loop 


0.4 


14000 


HB21 


0.8 


19000 


Gl 14.3+0.3 


0.7 


41000 


CTA1 


1.4 


24500 


HB9 


1.0 


7700 


Vela 


0.3 


11000 


G299.2-2.9 


0.5 


5000 


Monogem 1 


0.3 


86000 


Geminga 2 


0.4 


340000 



1 Plucinsky et al. 1996; 

2 caraveo et al. 1996; 



be assumed to be very much less than their distances from us. So, 
the models going to be adopted are the following : 

(1) For nearby sources (r<1.5kpc) : discrete both in space 
and time 

(2) For distant sources (r^l.5kpc): continuous both in space 
and time 



4 CR PROTONS IN THE SOLAR VICINITY 

Fig. 5 shows the model CR proton flux in the solar vicinity together 
with the observed data. The model spectra comprises of two com- 
ponents. One component (the distant component) is contributed by 
distant sources located beyond a distance of ~ 1.5kpc from us and 
the other component (the local component) is from nearby known 
SNRs located within a distance of 1.5kpc. The two components are 
treated seperately as discussed in the last section. 

From a compilation of several parameters of 231 all-type 
SNRs in the Galaxy based on radio observations as given in Table 
1 of Jian-Wen et al. (2005), it is found that there are only 1 1 SNRs 
which are located within a distance of 1.5kpc from the Sun. Out 
of these, 9 SNRs are with known ages. The two SNRs whose ages 
are not known are G31 1.5-0.3 and Lupus Loop, each located at a 
distance of ~ 1.2kpc from us. In addition to these SNRs detected 
in radio, there are also other SNRs within the region which are de- 
tected in other wavelength bands like Geminga and Monogem. The 
Geminga supernova is found to emit only X-rays and gamma-rays 
but not radio waves. Monogem is visible only in X-rays. The dis- 
tances from the Sun (r) and the ages (t) of the nearby SNRs with 
distances r < 1.5kpc whose ages are known are listed down in Ta- 
ble 1. The SNR RXJ1713.7-3946 is not included in the list because 
its distance is still not exactly known. The distance to the SNR is 
(6 ± l)kpc in contrast to the distance of lkpc estimated from soft 
X-ray obsorption (Koyama et al. 1997). Its age is approximately 
10 4 j/rs for a distance of r — 6kpc or 2 x 10 3 yrs for r = lkpc. 

The calculation of the distant component is done for a SNE 
rate of two times the Galactic value. It can be seen from the figure 
that non of the observed nearby SNRs nor the distant sources are 
able to contribute the observed data effectively below ~ 100GeV\ 
It is because for low energy protons, their diffusion time from dis- 
tant sources becomes comparable to their ionization, coulomb and 
adiabatic energy loss timescales due to which their fluxes are sig- 
nificantly suppressed. Also they may be convected away from the 
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Figure 3. Temporal density variations of 10GeV CR proton at a given 
location in the Galactic disk due to a source located at distances r = 
(0.1, 0.5, 1, 1.5)kpc (top to bottom) from the point. The calculation is for 
a particle release time of t = from the source. 
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Figure 4. Same as in Fig. 3 but for a proton energy of lOTeV . 



Figure 5. The flux of CR protons near the Sun produced by distant sources 
and nearby known SNRs listed in Table 1. The spectrum marked as G is 
the contribution from distant sources located beyond l.bkpc from us and 
the spectra numbered as (1,2, 3, 4, 5, 6, 7, 8) are the contributions from 
individual SNRs with 1— Geminga, 2— Monogem, 3— Vela, 4— Cygnus 
Loop, 5- Gl 14.3+0.3, 6- G299.2-2.9, 7- HB21 and 8- G65.3+5.7. The 
calculations are done for total output source energy W = 10 51 ergs, T = 
2.45, a = 0.6, D(WGeV) = 10 2S cm? s' 1 and V = l^kms^kpc' 1 
assuming a burstlike injection of the protons at time t = 0. The Solar modu- 
lation parameter <E> = 400M V . Data points shown corresepond to the local 
fluxes of CR protons measured by different groups (see Thoudam 2006 for 
details). 



Galactic plane before reaching us due to the presence of Galactic 
wind. Only high energy protons which diffuse faster can reach us 
from distant sources before they are lost in the intervening ISM. On 
the other hand, since the fluxes from nearby SNRs depend strongly 
on their distances r as well as their ages t, the low energy pro- 
tons from the knwon nearby sources have not yet reached us ef- 
fectively except those from Monogem and Geminga. But, below 
WOGeV the contributions of these two SNRs are approximately 
one to two orders of magnitude less than the observed data. Above 
WOGeV upto the maximum energy considered here (lOOTeV^), 
most of the contributions are from the distant sources. Among the 
nearby SNRs, only the Vela and Cygnus Loop are found to con- 
tribute the maximum flux in the high energy region. However, their 
contributions are around an order of magnitude less than the ob- 
served flux. 

To explain the observed proton spectrum below WOGeV, 
the presence of atleast one old SNR with an approximate age of 
t ~ 10 5 yrs located at a distance of r ~ O.lkpc is essential. The 
basic parameters (r, t) of the missing SNR is possible to be deter- 
mined if the particle spectra at different times (as shown in Fig.l) 
are carefully examined together with the variation of the maximum 
flux with r (shown in Fig.2). Since the maximum flux depends 
strongly on the distance of the source r from the Sun, the miss- 
ing source cannot be located far beyond ~ O.lkpc or nearer than it 
otherwise the proton flux will be either below or above the observed 
flux. Also, since the proton spectrum depends strongly on time t, 
the missing SNR cannot be younger nor older than ~ 10 yrs. In 
short, there is only one set of possible combination of (r, t) for the 
missing source which can explain the observed the proton spectrum 
below lOOGeV. 

Fig. 6 shows the predicted total local CR proton spectrum 
which consists of the contributions from distant sources located 
beyond 1.5kpc and that from an SNR with a characteristic age of 
t — 10°yrs positioned at a distance of r = 0.09fcpc from the Sun. 
The total spectrum does not include the contributions from the ob- 



Cosmic-Ray proton spectrum below 100 TeV in the local region 



— " 0.001 




local CR electrons where a nearby SNR is required to explain the 
observed high energy fluxes of electrons. 

Regarding the anisotropy that will be observed on Earth due 
to the missing source, the degree of anisotropy can be calculated 
using the relation 



Kinetic energy E (GeV) 

Figure 6. Total CR proton spectrum in the Solar vicinity. The thick solid 
line represents the total model spectrum which consists of the contribution 
from distant sources located beyond 1.5kpc from us (marked as G) and the 
contribution from a missing source (probably SNR) having a characteristic 
age of t = 10 5 yrs located at a distance of r = 0.09fcpc (marked as M). 
The calculation is done assuming burst-like injection of CR protons from 
the source at time t = 0. All the model parameters and the proton data 
points are the same as in Fig. 5 



served SNRs listed in Table 1 since their fluxes are found to be 
almost one-two orders of magnitude less than the observed flux. It 
can be seen that now the observed flux is explained effectively with 
the inclusion of the missing source. 



5 RESULTS AND DISCUSSIONS 

The model presented in the paper involves the seperate treatment 
of the distant and nearby sources for CR protons near the Sun. The 
distant sources are treated in a model continuous both in space and 
time whereas the nearby sources are discussed in a discrete model 
in space as well as time. The study shows that the total contribu- 
tion of the distant sources together with that of the observed nearby 
SNRs listed in Table 1 is not able to explain the observed CR pro- 
ton flux below ~ WOGeV. It is unlikely that it can be explained 
by the inclusion of the SNRs G31 1.5-0.3 and Lupus Loop (whose 
ages are not known) in the analysis since the CR proton flux de- 
pends strongly on the source distance r and both of them are at a 
distance of r ~ 1 .2kpc away from us. From fig. 2 it can be seen that 
the maximum proton flux from these two sources at WGeV ener- 
gies will be approximately 3—4 orders of magnitude less than the 
observed flux. To explain the observed CR proton spectrum below 
~ WOGeV, the presence of at least one nearby source with a char- 
acteristic age of t ~ 10°yrs located at a distance of r ~ O.lkpc 
is essential under the assumption of a burst-like injection of parti- 
cles at time t — 0. This missing source may be any old SNR which 
is unable to be detected by the present radio/x-ray telescopes be- 
cause of their limited sensitivities. In fact, assuming adiabatic phase 
the surface brightness estimated from an SNR of age ~ 10 yrs is 
fainter than the detection limit of the radio telescopes in studies of 
Galactic SNRs distribution (Kodaira 1974, Leahy & Xinji 1989). 
Such old SNRs can contribute a significant flux of CR protons in 
the solar vicinity although the electrons present inside it are not 
energetic enough to produce a strong radio/x-ray flux that can be 
detected by the present day telescopes. Incidentally, the parameters 
(r, t) of the missing SNR obtained in this work are found to be al- 
most same with that obtained by Atoyan et al. 1995 in the study of 



3 n 
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(Shen & Mao 1971) where r» and ti are the distance and age of 
the source respectively. For the source parameter obtained here i.e. 
Ti ~ O.lkpc and ti ~ 10 5 yrs, the degree of anisotropy of CR 
protons on Earth is found to be ^ ~ 5 x 10~ 3 on the assumption 
of burst-like injection of protons at time t — 0. 
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